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SUMMARY . yj 

A detailed mathematical proof is given in this report for the following 
new result: In the Restricted Three Body Problem with small mass ratio there 

exist one-parametric analytic families of synodically closed solution curves, 
which are near rotating Keplerian ellipses with rational sidereal frequencies 
and appropriate positive eccentricities. 

I. INTRODUCTION 

The equations of motion for the plane Restricted Three Body Problem can be 
written in the form 

x"+2ix'-x = -(1-p) (x-+p) |x4p | _3 - p(x4p-l) |xdp-l|“ 3 , (' = d/dt) (1) 

where x = x!+ix 2 is the complex position vector of the infinitesimal body P refer- 
red to a co-system rotating with angular velocity 1 about the center of gravity S 

of the two attracting bodies E and M of masses 1-p and p(0— p — 1) ns origin. 

When p = 0, the solutions of (1) are well known and can be represented as 
x(t) = e -it z(t), where the complex position vector z(t) describes the Keplerian 
motion: that is a solution of z" = -z|z|' 3 . Under suitable initial conditions 

this latter motion will be periodic. For instance, with 

z (0) = a(l+e) , z ' (0) = ic*/z(0), c* 2 = a(l-e 2 ), (s>0, 0<e<l ) , (2) 

z(t) moves along an ellipse with major half axis a and eccentricity e, having 
z = 0 as focus and z (0) at maximum distance from 0. Its sidereal period is 


To = 2jt | a 3 / 2 | . The corresponding x(t) will be periodic, iff. T 0 is commensu- 
rable with 2jt, or a “/ 2 = m/k, where k and m are relatively prime integers, 
m >0 and k is chosen positive respectively negative, if z(t) is direct respec- 

Vc 

tively retrograde, that is, sign k = sign c . The synodical period on the 

rotating ellipse then is T = 2jtm and the curve x = x(t), (0< t < T*) is closed 

after k-m positive revolutions around the origin. We denote this solution of 

(1) with p = 0 from now on by x (t) and obtain from (2) for its initial values 

* * . * *2 * 

X = a(l-fe), dx /dt = i (c -x ) /x at t = 0. (3) 

We shall show the existence of periodic solutions x(t) of (1) for small 

p > 0, which are near the generating solutions x (t) belonging to arbitrary 

integers k,m +0 and properly restricted e ; Namely, there are for fixed a = 

(m/k) 2 / 3 at most finitely many e in 0 < e < 1 with e = f (l-a“ 3 ) or with 
-k 

X (t) = 1 at least once in 0 ^ t ^ T . For every closed e- interval I, contain- 
ing none of these exceptional values, there exists a positive p* such that (1) 
possesses for every fixed p in 0 < p < p* <; 1 a family of periodic solutions 
depending analytically upon the parameter e in I. These solutions are holomor- 

phic also in p and transfer into x (t) for p = 0. Their synodical periods and 

Jacobi constants are holomorphic in e and p and depend both actually upon p. 

This result includes especially the existence of the periodic solutions 
of the so-called second kind for the Restricted Three Body Problem. Their 
existence had been claimed with supposed proofs by H. Poincare [5], K. 
Schwarzschi Id (jj] and C. L. Charlier [j2 J, whose invalidity was shown by P. 
Staeckel |YJ and A. Wintner [9], [lo], however. In these attempts the con- 

tinuation method of Poincare was employed in an isoperiodic or an isoenergetic 
manner. We too apply this continuation method, but replace the general per- 
iodicity condition of Poincare with a more special one which is based on the 


2 



symmetry of the dynamical problem (1) , and has already been used by G. D. 

Birkhoff [l ]. This condition is not only simpler and more natural, it also 
points out and reduces the redundancy in the classical periodicity conditions, 
which caused the critical functional determinants to vanish (Ref. ^9j and [?])• 
Otherwise, we achieve our goal by employing appropriate variables, which render 
the dependence of the Keplerian motion upon its initial values in a most simple 
form. 

In this regard, it is to be mentioned that G. D. Birkhoff j^L J showed, 
among others, the existence of periodic solutions of (1) for small p, which 
close in the rotating co-system only after sufficiently many revolutions about 
the mass 1-p, if their Jacobi constant determines a simple closed zero-velocity 
curve containing the orbits in its interior. B. 0. Koopmann j^3j established 
the analogon for the exterior case of a zero-velocity curve with forbidden 
bounded interior. More recently, J. Moser [4] showed the existence of periodic 
solutions of (1) for small p, which close after many revolutions and are near 
solutions of the existing first kind of Poincare (generated from circular 
motions for p = 0) . All these solutions correspond for p = 0 to periodic 
motions along rotating ellipses with rational a 3 / 2 . Presently, it is not 
known if these solutions for p > 0 coincide with certain of our above solutions 
x ( t) generated from x (t) with large |k-m| and suitable, small e > 0. 

Finally, we remark that several of these solutions for different m/k, whose 
existence is shown here for small p > 0, have been numerically calculated by 
us for increasing p on high speed electronic computers. They are particularly 
of interest when a(l-e)< 5 and 1 < a(l+e) < 1+& with small 5 > 0, since then 
they pass repeatedly near both masses of the Restricted Three Body Problem. 

The calculations indicate their existence for values of p at least as large as 

3 


that for the case 1-p/p. equal to mass of the Earth/mass of the Moon. Thus, 
thsir practical significance for astronautics is apparent and was one of the 
incentives for the investigation presented here. 

II. EXISTENCE PROOF 

Let x(t) be a solution of (1), which is holomorphic on an interval 
0<t^T, T > 0; that is, free of collisions. If (with a bar denoting the 
conjugate complex number) 

x(i T ) = x(^T) , x'(iT) = -x'(iT), (4) 

then the function x(T-t) of t is identical with x(t)^ since it also satisfies 
(1) and the two functions and their first derivatives coincide respectively at 
t = ^T. This implies^ if additionally 

x(0) = x(0) , x 1 (0) = “X ! (0) y ( 5 ) 

that x(T) = x(0), x f (T) = x 1 (0 ) } so that x(t) will be periodic with period T, 

since (1) is autonomous. Then the closed curve x = x(t), (0 < t < T) is sym- 

metric over the x^-axis^ since x(-t) = x(t). Especially x (t) satisfies (5) 
by (3) , and also (4) with T = T =|k|T Q , since as a consequence of (2) z(k-T Q /2) 

is real and z'(k»T 0 /2) is pure imaginary for every integer k. 

We introduce new real variables F, H, U, V instead of x = x x + ix 2 and 
y = x' + ix = y x + iy 2 by 

F = arctg X 2 / Xl , H = i(yf + y|) - r " 1 - c, r = (xf + x§)^, 

( 6 ) 

l u = x x /r - cy 2 , V = x 2 /r + cy x , c = x x y 2 - x 2yi . 

The functional determinant of this transformation is D = - (2H+3c) • c 2 /r 2 
Now c and H are first integrals of (1) for p = 0, since with z(t) = e ll: x(t) 
then 
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c = Im xy = Im zz 1 = const, of area, 

H + c = 'i|z T | 2 - | z j “ 1 = const, of energy, 
as well known for the Kepler Ian motion. By (2) then 

c = c , H = -c - 1 /2a , r = | x (t) |;> a(l-e)> 0 on x"(t). 

Thus 

Vc ic 2 

D = -(c - 1 /a) c /r 2 on x (t) , (7) 

and the transformation (6) is analytic and locally 1-to-l in a neighborhood of 

every point on the trajectory x (t), y (t) = dx (t)/dt + ix"(t), (0 < t < T ) 

* 

if D =f 0 or ac f 1, This holds for 0 < e < 1 always, when a < 1 or when 
e 4 v r (l-a- 3 ) for a > 1, c > 0, and this assumption will be made from here on. 
Now (6) transforms (1) in case \i = 0 into 
F 1 + 1 = c/r 2 = c“ 3 (1-UcosF-VsinF) 2 , 

* ( (8) 

H' = 0, U 1 = V, V' = -U, 

since from (6) c = r (1-UcosF-VsinF) /c. If we denote initial values by the 
corresponding small letter, (8) can be integrated up to 
H = h, U +iV = (u+iv) e 

' , (9) 

F + 1 = c“ 3 (l-ucos (F+t) - vs in (F+t) ) 2 . 

This and (6) yield 

u 2 + v 2 = U 2 + V 2 = 1- 2c 2 /r + c 2 (yf + y%) = 1 + 2hc 2 + 2c 3 , (10) 

so that c depends upon the initial values h, u, v. Substituting for these the 
special values 

+ ★ it -k 

f = 0, h = -c -l/2a , u = €, v =0 (11) 

•>V 

derived from (3) and (6) gives c = c and by (9) the original solution x (t), 
but now represented in the new variables. 


5 



In general (6) transforms (1) into a system 

F' = gl , H 1 = g 2 , U f = g 3 , V* = g 4 , (12) 

where the g^ = g n (F,H, U, V, p) are holomorphic functions of all variables in a 

neighborhood of the special solution determined by (9), (11) and for suffi- 
ce 

ciently small p > 0, if x (t) is free of collisions for \i - 0. Given a = 
(m/k) 2 / 3 , such a collision, that is, x (t) = 1 for some t in 0 < t < T , can hap- 
pen only for finitely many values of e in 0 < € < 1. We return to this con- 
dition later and assume here only that these e are omitted. Then, according 
to Poincare’s extension of Cauchy's existence theorem for ordinary differential 
equations (for a modern proof see Ref. [7]) the solutions of (12) are holomor- 

k 

phic functions of t,f,h,u,v, and p for 0 < t < 2T , say, and sufficiently small 

| f | + |h-h*| + |u-e| + | v | + (jl, p > 0, using (11). For p = 0 (12) becomes (8). 

We consider now the solutions of (12) with small p > 0 and initial values 

near (11). These solutions can be assumed remaining near the original solution 

* 

given by (9), (11) for 0 < t < 2T , so that especially (6) is applicable. In 
order that such a solution will be periodic in the former coordinates x x , x 2 
with period T > 0, it is sufficient according to (4), (5), and (6) that 
[f C^T, f , h, u, v, |i) = jt(k-m) , f = 0 , 
lv(£r,f,h,u,v,n) =0 , v = 0 , 


(13) 


since for the considered solutions of (12) c remains near c , and therefore 

1 * 

the assumption c f 0 at t = 2 T can be made, if T < 2T . These equations are 

k 

actually satisfied for the original solution x (t) and p = 0, or with (11) for 
k k k k k 

T = T = 2mjt, f=f, h = h , u = u, v = v , |i = 0 , (14) 


since 


k i k 

arc x ( 2 T ) = arc 


fe im7t z ( j k | T q / 2) j = arc|e J _ (k-m) n 
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by (2). Hence (13) can be satisfied for small p > 0 and initial values near 
(14), if for instance the functional determinant 

D* = F t V h - F h V t +0 for t = |l at (14). (15) 

By the holomorphy it suffices putting p = 0 in F,V first and then calculate 
the partial derivatives, which therefore can be found from (9) and (10). De- 
noting partial derivatives by a corresponding index, (9) gives = 0 and 

V = -U, hence with (11) and (14) 

■fr k 

D = e cos mn • F h (mjt,0,h , e,0), (p=0). (16) 


We put 

* = t + F(t,0,h, e,0) = 4>(t,h) . (17) 

Then it follows from (9) that 4> is uniquely determined by inversion of the 
integral 


t = c 3 / (1-e cos\[r)~ 2 d\|r . 


(18) 


The value of c here is determined by e 2 = l+2hc 2 +2c 3 from (10). This implies 
c h = -c (2h+3c) " 1 . (19) 

Differentiating (18) partially with respect to h gives 

0 = 3c 2 c^ I (1-e cos)|0 ~ 2 d\]/ + c 3 (l-e cos d)” 2 4>^ , 


hence with (19) * 

4 > 

k k k k P 

4*. (mjt, h ) = 3 (1-e cos 4 1 ) 2 (2h 4- 3c ) -1 / (1-e cos \|/) _2 d\jr , 


( 20 ) 


k k 

where 0 - 4>(mjr, h ). But 


*3 


2c / (1-e cos \|r)" 2 d^ = +T 0 = 2jtm/k 


for the period on the original Keplerian ellipse as well known, thus by (18) 
4> = kjt and by (16), (17), and (20) finally 
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( 21 ) 


D* = 3e(-l) m (l-e(-l) k )%it/ (c^a-^c* 3 * 0 , 
so that (15) actually holds, since 0 < € < 1 and c =t= a" 1 as required for (7). 
From the implicit function theorem follows now that (13) can be solved 

JU JU 

for T and h in a neighborhood of (14), and that T-T' and h-h for u=u =e result 
as power series in p without constant terms having positive radius of con- 
vergence. This implies especially 0 < T < 2T for small p > 0 as assumed be- 
fore, so that by (13) the existence of the desired periodic solutions x(t) is 
now actually shown. Their initial values x(0), x 1 (0) are determined by (6) and 
f = 0, h = h(p, e) , u = e, v = 0 (22) 

as functions of p, e and a. 

Since for any given a = (m/k) 2 / 3 the foregoing is valid as long as € is 

not one of the previously excepted values, it is a consequence of the local 

existence theorem for implicit functions and of the covering theorem that 

T = T(p, e) and h = h(p, e) are holomorphic functions of p and e on every closed 

e-interval not containing one of the exceptional values and on 0 ^ n < n 

ic 

with corresponding sufficiently small p > 0, The solutions of (12) belonging 
to (22) are holomorphic functions of their initial values and thus the corres- 
ponding periodic solutions x(t) are holomorphic on 0 < t < T and in e and p as 
stated in the introduction. 

If one introduces in the transformation (6) instead of the variable H the 
Jacobi integral 

J = !|y | 2 - c -(i-i-O I x+ m-I~ 1- ^l x 4 u _ 1 l _1 ? 

our whole consideration can be carried through in the same way. Thus also 
J = J(p,e) is a holomorphic function of both variables, which clearly follows 
from 
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J = h + l x 0 l _1 “ ( L-m.) | | " 1 - |i|x o -+^i-l| 1 , X q = x (0) 

and the foregoing, too. In (13) then h is to be replaced by the initial value 
j of J, but J = j now. It is then of interest to consider besides (15) the 
other two functional determinants suggested by (13): namely, 

-j- 

Dt = F V - F V and Do = F .V - F V . 

± tu Ut JU U J 

1 “ 

for t = ^T at (14), where j = h . Since here = -sin mic = 0, V = 0 and 

* * Jc 

V = -U by (9), it follows D 2 = 0 and D x = e cosmje (mjt, 0, h , e,0). 

Putting 

^ = t + F(t,0,h ,u,0) = vlr(t,u) , 

it follows from (9) \|r = c _3 (l-u cos \|/) 2 and thus similar as from (18) for $ 

now \|r(mjt, u) = kjt identically in 0 < u < 1. Hence also D x ~ 0. Thus, solvabil- 
ity of (13) with respect to T, u for fixed j = h , or with respect to j,u for 

■JL 

fixed T = T remains at least doubtful, if at all possibla In fact, A. Wintner 
( 9^ has shown that for sufficiently small e > 0 isoperiodic solutions do not 
exist. 

It is interesting to note that the treatment presented here is also 

effective for the periodic solutions of the first kind. In this case, e = 0 

and a 3 / 2 = + T Q /2jt = w” 1 can be taken arbitrarily, especially not rational > 

with the exceptions u) 4 -1,1 and w 4 1 + rrT - 1 , m natural. Then it follows 
* * * 

Di 4 0; D 2 ^ 0, but D = 0 in (15), which is the reverse of our present situ- 
ation for e > 0. Of course, an existence proof for the solutions of the first 
kind is long known. 

Finally, we consider the restrictions placed upon e for our existence 
proof. These are 
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( 23 ) 


0 < e < 1, e 4 /(1-a -3 ) for a > 1, c* > 0, 

x (t) 4 1 in 0 ^ t ^ mjt for a = (m/k) 2 / 3 . 

The first of them is equivalent to ac f 1 or 2h + 3c 4 0 and was required 
for (7) and (21). The dynamical meaning of this condition becomes clear, when 
the sidereal frequency w = 2jt/T Q = |a" 3 / 2 | is used instead of a to characterize 
the generating elliptic motion together with its eccentricity. Since here 
p = 0, H coincides with the Jacobi- integral J, and then (19) shows that the 

k 

Jacobi-constant j = h of the generating motion x (t), for which now 

j = -c - l/2a = - •ga) 2 ^ 3 - sign c » w -1 / 3 /(1-e 2 ) , (w > 0) 

, * 

has for given e as function of c or w a relative extremum (absolute maximum) 
k „ 

for ac = 1 or u - /(1-e 2 ), in the direct case only. The inverse function and 
thus the sidereal period T^ then cannot be defined uniquely as functions of the 
Jacobi constant and the eccentricity in a neighborhood of the branch points 
(j* = -(1-e 2 ) 2 / 3 • 3/2, e). The dynamical meaning of (23) is to exclude collis- 
ions with the perturbing body. 

We shall show that for fixed k and m (23) excludes at most finitely many 

values of e. To see this, we represent the Keplerian motion z = z(t) = 
it * 

e x (t) with the help of the eccentric anomaly w, namely with (2) in the 
well known form 


z = a(e + cos w + i/(l-e 2 ) sin w) , 
t - a 3 / 2 (w + e sin w) 


If (23) does not hold, then with appropriate e between 0 and 1 

z = e lt: , |z| = a (1 + e cosw) = 1 (23) 

for some w between 0 and k n inclusively, using (24). The last equation implies 
1 1- 1/a |^: € < 1, thus a > j|. Hence (23) holds always, when 2a< 1. If a = 1, 
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then by (25) cos w = 0 and thus with (24) 
e = cos t = cos(wj^) = - sin 

which is impossible for 0 < e < 1. Let now a > -|, a f 1 and (25) be satisfied 
for appropriate e and w. Then with -A = 1-1/a and (24) 
e cos w = A £ 0, a(e + cos w) = cost, 

a(cos w 4- A/cos w) - cos|a 3 ^ 2 (w + Atg w) = 0. (26) 

Since ] A j < e < 1 and with w* from cos w* = |a|, 0 < vrkjt/2, we have here 

|w - njt | < w' c , n = integer, |n|<; |k|. (27) 

Now the left side in (26) is a holomorphic function of w in (27) and even in 
all of the f ini t ely-many closed circles |w-njt|< w* . Thus it has at most 
f ini tely-many zeros in (27), and these correspond to at most f ini tely-many 
e = A/cos w satisfying (25). This proves our statement about (23). 

One can easily see that our whole derivation remains valid, if we begin 
in (2) with z(0) = a(l-e) at minimum distance from 0. This merely replaces e 
by -e in the subsequent equations. But it leads to new periodic solutions 
of (1) for small p, which are actually different from the previous ones, even 
when considering solution curves only, if k-m is odd. Additional such solutions 
are found starting with x~ = -a(l-fe) in (2), if k-m is even. All of our peri- 
odic solution curves are symmetric over the x-^axis, and they can be readily 
visualized in the rotating co-system. 

Added are a few illustrations calculated for the case p = 1/82: Figure 1 

and Figure 3 show synodically closed solution curves of (1) for m=l, k=2 and 
m=2 , k=5 respectively, in the rotating co-system. Figure 2 shows for the first 
case the circular paths of E and M and the path of P in an inertial co-system 
with origin at S. Figure 4 shows for the second case the paths of P and M in 
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a co-system with origin at E and space fixed orientation. Note the short dura- 
tion capture of P by M with subsequent rejection, which happens after elapse 
of every T units of time; i.e., every time when P has completed nearly k 
Keplerian elliptic orbits with focus at E. 
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FIGURE 2. CIRCULAR PATHS OF E , M AND THE PATH FROM FIGURE 1. IN 
INERTIAL CO. -SYSTEM FROM t=0 TO NEARLY t=2T. CAPTURE 
AND REJECTION OF P BY M. 
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FIGURE 
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FIGURE 4. CIRCULAR ARC OF M, AND THE PATH FROM FIGURE 3 AS SEEN 
FROM E UNDER SPACE FIXED ORIENTATION NEARLY FROM 
t=0 TO t=6/5 T. CAPTURE AND REJECTION OF P BY M. 
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